Ab initio study of the migration of intrinsic defects in 3C-SiC 
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The diffusion of intrinsic defects in 3C-SiC is studied using an ab initio method based on density 
functional theory. The vacancies are shown to migrate on their own sublattice. The carbon split- 
interstitials and the two relevant silicon interstitials, namely the tetrahedrally carbon-coordinated 
interstitial and the (llO)-oriented split-interstitial, are found to be by far more mobile than the 
vacancies. The metastability of the silicon vacancy, which transforms into a vacancy-antisite complex 
in p-type and compensated material, kinetically suppresses its contribution to diffusion processes. 
The role of interstitials and vacancies in the self-diffusion is analyzed. Consequences for the dopant 
diffusion are qualitatively discussed. Our analysis emphasizes the relevance of mechanisms based 
on silicon and carbon interstitials. 

PACS numbers: 61.72. Ji,66.30.Lw,66. 30. Hs 
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I. INTRODUCTION 

Silicon carbide is an important wide band gap semi- 
conductor which offers a variety of applications for high 
power and high frequency devices. During device pro- 
cessing, various intrinsic defects are introduced that af- 
fect the electronic properties of the material. The mobile 
intrinsic defects, namely vacancies and interstitials, play 
a pivotal role in self diffusion and dopant diffusion as well 
as in the annealing of ion-implanted material. 

Substitutional impurities are per se immobile. They 
need vacancies or interstitials as vehicles for the migra- 
tion. The rapid diffusion in the course of annealing can 
considerably affect the implanted dopant profiles. For 
example, a transient enhanced dopant diffusion was ob- 
served in boroni£ and aluminum 3 implanted 4H-SiC. It 
is initiated by an excess concentration of intrinsic de- 
fects, most likely silicon interstitials as the analysis of re- 
cent experiments indicatesiSaiii in contrast to the earlier 
assumption 6 of a vacancy-mediated mechanism. Similar 
mechanisms are operative in the more fundamental self 
diffusion. As a matter of fact, a microscopic understand- 
ing of the mechanisms underlying the self diffusion should 
provide insight into the mechanisms of the dopant diffu- 
sion. In the most recent studiesi2*2ii& of self diffusion in 
3C- and 4H-SiC diffusion constants have been measured 
for the carbon and silicon self diffusion. The basic diffu- 
sion mechanisms, however, have not been unraveled. In 
particular, the role of interstitials and vacancies in differ- 
ent diffusion processes remains open for both polytypes. 

Alongside with the diffusion, the annealing of mo- 
bile intrinsic defects contains information about the mi- 
gration of vacancies and interstitials. The annealing 
of vacancy-related defects has been studied by positron 
annihilation spectroscopy (PAS)ii*i2ii 3 ji4 and electron 
paramagnetic resonance techniques (EPR)A£4i£iiLi& in ir- 
radiated material. The identification of the EPR-centers 
as isolated siliconi£ii§i and carbonic vacancies has been 
verified theoreticallyiiSiiS It was found that the anneal- 
ing behavior of carbon and silicon vacancies shows strik- 
ing differences. For example, in EPR-Experimentsi^ the 



silicon vacancy was found to anneal in several stages at 
150°C, 350°C and 750°C, whereas for the carbon vacancy 
a single characteristic annealing temperature of 500° C 
was observed^ Recently, an EPR-center was observed 
by EPR experiments^ in material annealed above 750° C 
and interpreted as a vacancy-antisite complex. This com- 
plex may be a decay product of the metastable silicon va- 
cancy (in p-type and intrinsic material) as predicted by 
theory^i^ A comprehensive interpretation of annealing 
experiments in terms of elementary processes crucially 
depends on an understanding of the underlying diffusion 
mechanisms of intrinsic defects. 

Theoretical investigations by ab initio methods can 
provide important insight in a microscopic picture of 
the self- and dopant-diffusion as well as the annealing 
kinetics. In the present paper we investigate the mo- 
bile intrinsic defects and their migration mechanisms in 
3C-SiC using a method based on density functional the- 
ory. The Frenkel pair recombination and other diffusion- 
controlled annealing mechanisms are treated elsewhere*^ 
Our investigation is based on recent theoretical stud- 
ies of the ground state properties of vacancies^S^ 
and antisites^ and addresses interstitials and vacancy- 
antisitc-complcxcs that are relevant for the migration. 
We show that different types of silicon interstitials, which 
have entirely different migration paths, are important in 
p-type and n-type material. The metastability of the sil- 
icon vacancy in p-type material is demonstrated to have 
strong implications for vacancy-assisted diffusion mech- 
anisms. The considered diffusion processes are summa- 
rized in Fig. ^ which shows the most important migra- 
tion channels for vacancies and interstitials in 3C-SiC. We 
find that interstitials are by far more mobile than vacan- 
cies and therefore play a prominent role in diffusion. We 
also demonstrate that the diffusion mechanisms strongly 
depend on the charge state of the mobile defects. 

The outline of the paper is as follows. Beginning with 
a description of the method in Section ^ we turn to 
the ground state properties of interstitials, vacancies and 
vacancy-antisite-complexes in Section IIIII and IIVI The 
migration mechanisms of vacancies and interstitials are 
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FIG. 1: The dominant migration mechanisms of interstitials 
and vacancies: (a) carbon vacancies (Vc) migrate on the car- 
bon sublattice only, (b) transformation of the silicon vacancy 
(Vsi) into a carbon vacancy-antisite complex Vc-Csi (p-type 
and compensated) and the vacancy migration on the silicon 
sublattice (n-type), (c) carbon interstitial migration via split- 
interstitial configurations C sp (ioo) and C ap si(ioo> , an d (d) kick- 
out mechanism for the carbon coordinated silicon interstitial 
(SiTc) via split-interstitials (Si sp (i o> and Si sp (n ))(p-type) 
and direct migration of the silicon split-interstitial Si ap (no) 
(compensated and n-type material). 
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FIG. 2: Madelung corrections for the silicon interstitial Si^c- 
The formation energy i5si,f as calculated in the 64 and 216 
atom cell with and without the Madelung correction (filled 
and open symbols, respectively) is plotted versus the inverse 
defect-defect distance. The squares and circles refer to a 
(2x2x2-mesh ) k-point set and the F-point. The extrapo- 
lation of £!^ lTC obtained in the 216 atom cell to the isolated 
defect using the monopole correction is shown by the solid 
line. 



analyzed in Section The role of interstitials and va- 
cancies in the self- and dopant-diffusion is discussed in 
Section IVII A summary concludes the paper. 



II. METHOD 

We carried out first principles calculations using the 
plane wave pseudopotential package FHI96MD 26 within 
the framework of density functional theory (DFT)i2L2£ 
The local density approximation 29 (LDA) is employed for 
the exchange-correlation functional and spin effects are 
included within the LSDA where necessary. 

We describe the intrinsic defects and their environment 
using periodic supercells. Large supercells with 64 and 
216 crystal lattice sites are used. The electron transi- 
tions between the defect and its periodic images give rise 
to the formation of defect bands and affect the defect 
energetics. This artificial defect-defect interaction can 
be efficiently reduced by a special k-point sampling^ A 
special k-point mesbJi with 8 k-points in the Brillouin 
zone (2x2x2-mesh) provides converged defect energies 
for the 64-atom cell. It is crucial to preserve the correct 
degeneracy of the isolated defect levels, as the degener- 
acy may lead to a Jahn- Teller instability. In those cases, 
where symmetry-lowering Jahn- Teller-distortions are im- 
portant, we perform the calculations using only the T- 
point in a 216 atom cell (for which already sufficiently 
converged formation energies are obtained), with the ex- 



act point symmetry of an isolated defect. In the case of 
charged defects we follow the approximate procedure of 
Makov and Payne^ 2 - to account for the electrostatic inter- 
action of the periodically arranged defects as well as their 
interaction with the compensating background. The rel- 
evance of such corrections has been already pointed out 
by Torpo et al. for the vacancies in 411-SiCj 2 ^ Especially 
for the highly charged interstitial defects this procedure 
improves the defect energetics considerably. In Fig. [21 we 
demonstrate this for the carbon-coordinated silicon in- 
terstitial Sixc (cf- Sec. lllljl . We find that the completely 
ionized interstitial silicon (charge state 4 + ) is strongly 
screened by the surrounding lattice. A comparison of 
the formation energies calculated for the 64 and 216 atom 
cell with and without the Madelung correction shows in- 
deed that the correction leads to a more consistent de- 
scription. Whereas we obtain 0.5 eV and 1.9 eV without 
the correction using special k-points, the corrected re- 
sults of 4.23 eV and 4.47 eV agree to within 0.24 eV. Here 
we have only included the dominant monopole correction 
(q 2 a/2eL, where a is the Madelung constant of the sim- 
ple cubic lattice, L the defect-defect distance and e the 
experimental dielectric constant), the quadrupole correc- 
tion (~ qQ/eL 3 , where Q is the quadrupole moment) 
should lead to a smaller correction as suggested by the 
detailed analysis of Lento et al£^ 

Soft norm-conserving pseudopotentials of the 
Troullier-Martins type^ and a plane-wave basis set 
are employed. The carbon pseudopotential has been op- 
timized for calculations with a small basis set (matching 
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radii rg=1.6Bohr, rp=1.7Bohr and rjj=1.5 Bohr). For 
the silicon potential the standard parameters are used. 
With these pseudopotentials and a basis set including 
plane waves of a kinetic energy up to 30 Ry defect 
formation energies have been calculated with a basis set 
induced convergence error below 0.1 eV (for an energy 
cut-off of 40 Ry all energy differences are converged) . For 
each defect the coordinates of all atoms in the supercell 
have been relaxed. This is essential to obtain a correct 
energetic description. 

A. Formation energies 

As a guide through the hierarchy of intrinsic defects 
one often uses their relative abundance in thermody- 
namic equilibrium. At constant volume the concentra- 
tion of defects is determined by the free energy of defect 
formation Ff. 

CD = csexp (-^y) ' (1) 

where cs is the concentration of the sites that are open 
to the defect, /cb is the Boltzmann constant and T is the 
temperature. The free energy of formation is specified by 
the formation energy Ef and the formation entropy St of 
the defect: 

Ft=E { -TSi . (2) 

Usually the defect abundance is well described by the for- 
mation energy Ef alone. Even though the formation en- 
tropy may amount up to 10 fcs (values of about 8 kn have 
been recently reported for vacancies in Si 3 ^ and GaAs2&), 
its contribution can be relevant only at fairly high tem- 
peratures and for defects with comparable formation en- 
ergies. 

In the compound material SiC the equilibrium de- 
fect concentration depends on the chemical environment 
which is characterized by the carbon and silicon chemi- 
cal potentials, /ic and fj,gi- The chemical potentials may 
vary only within certain bounds to prevent the forma- 
tion of other more stable phases than SiC. In excess of 
carbon (C-rich conditions) such a phase can be graphite. 
In the other extreme (Si-rich conditions), this would be 
crystalline silicon. The formation energy of a defect is 
calculated a a 37 i 38 

Ef = En, cell -ncMc- "Si A*Si - n c fi F , (3) 

where Ed, cc \\ is the total energy of the defect supercell, 
nc and nsi & r e the number of carbon and silicon atoms in 
the supercell. In case of an ionized defect, n c is the num- 
ber of excess electrons transferred from the Fermi-level 
/ip to the localized defect levels (for positively charged 
defects n e is negative). In equilibrium, the chemical po- 
tentials of the environment fic and fxgi are related to the 
chemical potential of the crystal /isic by //gic = Mc + Msi- 



With this relation we eliminate nc m Eq. ®. We ex- 
press the remaining free parameter nsi = Msi + ^A 4 m 
terms of the chemical potential of crystalline silicon /ig ; 
and the difference to this value A/i. These substitutions 
in Eq. Q yield 

£"" c = £"p c - (?i Si - n c ) A/i - n c , (4) 

where E^' = E£° coll - n c MSiC - ("Si - nc) Msi is ex- 
pressed by the number of the involved atoms and chem- 
ical potentials, which are obtained by total energy cal- 
culations for the perfect crystals. As described above, 
neither fxc nor /tgi can exceed the values of the corre- 
sponding bulk phases (graphite or diamond and silicon). 
Therefore A/i may vary only between the negative heat 
of formation of the corresponding polytype —Hf, sic & n d 
(C-rich and Si- rich conditions, respectively). Our cal- 
culated value of the heat of formation amounts to 0.61 eV 
and depends only slightly on the polytype. This result is 
in good agreement with the value of 0.58 eV obtained by 
Zywietz et alS^ and slightly lower than the experimen- 
tal value of 0.72 eV. Here we have deliberately chosen 
diamond as the reference. This choice is comparable to 
the alternative choice of graphite due to the small energy 
difference between graphite and diamond. Note that our 
definition of the formation energy En is independent of 
this aspect. The value of En of Ref. |2J is retained by 
evaluating Ef for A/i = — 1/2 Hf. sic- For the calculation 
of charged defects it is convenient to state the value of 
the Fermi-level /ip relative to the valence band edge Ey. 
In this case the energy En is given by the formula 

= En, con - nc MSic - ("Si - nc) Msi ~ n cEv , (5) 

The value of Ey is not reliably accessible in the calcu- 
lation of the involved defects. Therefore we take the 
value Ey from a calculation for the perfect crystal and 
remove shifts in the average potential by aligning the 
bulk-spectrum features in the density of states of the de- 
fect cell and the defect-free cell. The alignment is con- 
sistent within 0.06eV over the whole energy range. By 
this procedure we also verified that the defect cells con- 
tain a sufficiently large bulk-like environment around the 
defect. 

The realization of a specific charge state of a deep de- 
fect for a given Fermi-level /iF requires that the formation 
energy of the defect in the corresponding charge state is 
lower than for the other possible charge states and that 
all occupied single particle levels lie below the conduction 
band edge. Even though the Kohn-Sham levels have no 
physical interpretation by the design of the theory (only 
the electron density and the total energy are a well de- 
fined quantities), they reproduce the experimental quasi- 
particle band structure quite well. Yet, the calculated 
Kohn-Sham band gaps are usually smaller than experi- 
mental values (in 3C-SiC and 4H-SiC we obtain 1.2 eV 
and 2.2 eV as compared to the experimental findings of 
2.39eV and 3.27 eV). In this paper we follow a common 
practice (cf. e.g. Ref.|25J) and use the experimental value 
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for the conduction band edge. With this approach it 
is not possible to unambiguously determine the charge 
states that are only stable in a region close to the con- 
duction band edge. 

The thermodynamic ionization level of a charged defect 
is given by the value of the Fermi-level at which the defect 
alters its charge state. It is obtained from the formation 
energy by 

e(q i \q l )=E^-E§ . (6) 

Here qi and q2 indicate the different charge states of the 
defect. Usually only one electron is transfered between 
the electron reservoir and the defect levels. A simulta- 
neous transfer of two electrons is unfavorable due to the 
cost of the effective electron-electron repulsion. How- 
ever, in some cases^iS the electron-electron repulsion 
can be compensated by a configurational relaxation aris- 
ing from a strong electron-phonon coupling. This effect, 
known as the negative- U effect, leads to an attractive 
effective electron-electron interaction and the ionization 
level s(q — l\q) appears below e(q\q+ 1). 



B. Defect migration 

In order to analyze possible migration paths of intersti- 
tials and vacancies we apply two different standard meth- 
ods. For the interstitial migration we use an implemen- 
tation 41 of the ridge method of Ionova and Carter. 42 In 
this method, a saddle point search is conducted for given 
initial and final configurations of the migration event, 
both configurations being slightly distorted towards each 
other. The two configurations are iterated such that they 
approach each other and the energy barrier along the 
line connecting the two configurations is minimized. By 
this procedure both configurations converge to the sad- 
dle point and the lowest migration energy barrier E m is 
obtained. The search is repeated for all relevant charge 
states and the obtained transition states are analyzed. 
This automatic search of the saddle point fails when de- 
fect levels of different symmetry cross along the migration 
path. In this situation, which applies to the migration of 
vacancies, we analyze the potential energy surface using 
the drag method. This approach has been previously 
employed in a study of the gallium vacancy migration 
in GaAsi^li For a relevant set of coordinates (reaction 
coordinates), the potential energy surface is calculated 
by constraining the coordinates of interest and minimiz- 
ing the total energy for all remaining coordinates. The 
consistency of our choice of the reaction coordinates is 
verified. This includes the smoothness of the potential 
energy surface as a function of the reaction coordinates. 
Details of the procedure are given in Sec. |VJ where the 
individual migration mechanisms are discussed. 



III. INTERSTITIALS 
A. Silicon interstitials 

SiC as a compound semiconductor possesses vari- 
ous interstitial configurations, some of them being dis- 
tinguished only by a carbon- or silicon-like environ- 
ment. Important interstitial configurations are depicted 
in Fig. |3| There are tetrahedral configurations with either 
four silicon or carbon neighbors and a hexagonal con- 
figuration. Besides these, split-interstitial configurations 
and bond-center configurations can exist. In a split con- 
figuration an interstitial silicon atom shares a site with 
a lattice atom, which could be either carbon or silicon. 
These dumbbell-like interstitials occur with orientations 
of the atom pair in the (100)- and the (HO)-directions. 
In a bond-center configuration the interstitial silicon is 
centered at the bond of two neighboring lattice atoms. 

We have considered all of the above configurations in 
our investigation. The hexagonal and bond-center in- 
terstitials as well as the split-interstitials on the carbon 
sublattice turned out to be unstable in all relevant charge 
states. The calculated formation energies for the stable 
interstitials in comparison with the formation energy of 
the silicon vacancy (c.f. Sec. lIVp are displayed in Fig-EJas 
functions of the Fermi- level fj,p. To make the relation to 
Eq. Q more transparent, we consider the case A/x = 0, 
i.e. Si-rich conditions, in Figs. and |3J Since in case of 
silicon interstitials nsi — «c = lj their formation energies 
rise by 0.61 eV when referring to C-rich conditions, while 
the formation energy of the silicon vacancy drops by the 
same amount. The charge states of the defects in the 
different ranges of hf are indicated. Except for the split- 
interstitial Si sp (no) all silicon interstitials are positively 
charged. 

The hierarchy of silicon interstitials strongly depends 
on the position of the Fermi-level. This is an immedi- 
ate consequence of the existence of the deep interstitial 
levels. In p-type material the most abundant is the tetra- 
hedrally carbon-coordinated interstitial Sixc followed by 
the split-interstitial Si sp /no). Around mid-gap, this or- 
dering changes: as the interstitial Sixc is completely ion- 
ized (carrying a charge of 4 + ) its formation energy rises 
more rapidly than that of the Si sp (no) interstitial, which 
is neutral for /ip around mid-gap. For Sixc we can only 
estimate the formation energy of the charge states 3 + 
(c.f. below). This also applies for the split-interstitial 
Si S p(ioo) ■ Therefore we have plotted the formation energy 
only up to the estimated position of the lowest ionization 
level. The tetrahedrally silicon-coordinated interstitial 
Sixsi and the split-interstitial with (fOO)-orientation have 
a higher formation energy than the other stable configu- 
rations. 

In the following we describe the electronic structure 
of silicon interstitials. The interstitial Sixc does not in- 
troduce deep levels into the band gap. The silicon atom 
loses its valence electrons and becomes completely ion- 
ized. However, the interaction with the carbon and sili- 
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FIG. 3: Silicon-interstitials: (a) formation energy for Si-rich 
conditions in comparison to the silicon vacancy (b) geome- 
try of investigated interstitial sites; the hexagonal site is not 
stable. 



con neighbors induces an electron flow to the silicon in- 
terstitial, resulting in an electron distribution that re- 
sembles that at a regular silicon lattice site. The s- 
and p-components of the electron density integrated over 
a sphere with a radius of 1.75 Bohr (half the nearest 
neighbor distance in SiC) around the Si-atom deviate 
from the values for a silicon site in a perfect crystal by 
less than 1%. This is also true for spheres of smaller 
radii. In Fig. 0] we have plotted the charge redistribution 



5 n = n Si 4+ — n-buik, where n- Si 4+ is the electron density of 
the supercell containing the interstitial and nbuik refers 
to the density of the perfect crystal. The polarization 
of the surrounding carbon-silicon bonds and the forma- 
tion of bonds with the silicon interstitial is clearly visible. 
Although we have not found any localized one-particle 
levels within the band gap, we do not expect the charge 
state 4 + to persist throughout the whole experimental 
band gap. Indeed, slightly above the experimental con- 
duction band energy we observe a localized state. Due 
to the large charge state correction the corresponding 
ionization level e(3 + |4 + ) lies substantially lower, in a re- 
gion somewhat above the experimental mid-gap. Thus 
the charge state 4 + should be stable at least up to the 
mid-gap position of the Fermi-level. Yet, the quantita- 
tive determination of the ionization level is not possible 




FIG. 4: Binding of Sitc: polarization Sn(r) of the electron 
density around the bare interstitial silicon ion at the SiTC-site 
obtained as a difference between the electron densities of the 
interstitial nsi TC;4+ and of the bulk crystal nbuik (calculated 
in the 216-atom cell). 



in DFT-LSDA calculations due to the well-known band 
gap problem. The localized level is in resonance with 
the Kohn-Sham (and experimental) conduction band. To 
obtain their relative position (which is crucial for a cor- 
rect occupation of levels) one would have to apply a XC- 
discontinuity correction to the conduction band and the 
localized state, which is well beyond the scope of the 
present work. This problem is less acute in 4H-SiC, where 
the band gap is substantially wider. We found^ 3 . that 
the carbon-coordinated Si-interstitial in 4H-SiC possesses 
a similar ionization level e(3 + |4 + ) and a corresponding 
one-particle state that is located within the Kohn-Sham 
band gap. 

The screening of the silicon interstitial weakens the 
neighboring carbon-silicon bonds. This is manifested by 
an inward relaxation of the carbon neighbor-shell to- 
wards the interstitial by 3.5% of an ideal bond length. 
The energy associated with this relaxation amounts to 
4.82 eV. The relaxation of the first and second neighbor 
shell yields 4.37 eV, whereas the inclusion of only the first 
neighbor shell results in a gain of 0.34 eV. This behavior 
shows the strong coupling between the induced polariza- 
tion and the relaxation. 

The tetrahcdral interstitial Sixsb m contrast to SixC) 
has a deep non-degenerate level in the band gap below 
mid-gap. Its occupation enables charge states between 
4 + and 2 + . However, due to the charge state correc- 
tions, Sirpg; is not realized. Only a single ionization level 
(2+|3+) at 0.2eV is present (c.f. Tab. |T|). The state has 
s-like character at the interstitial silicon and p-like char- 
acter at the fourth neighbor carbon-atoms, which have a 
common bond with four silicon neighbors of Sixsi- There 
is only a negligible amplitude at the first neighbor shell. 
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TABLE I: Ionization levels of carbon and silicon interstitials 
with respect to the valence band edge Ev in eV. The calcu- 
lations have been performed in a 64-atom cell using special 
k-points (C sp (ioo) and C sp si(ioo> : 216-atom cell and T-point to 
include spin-polarization and Jahn- Teller distortions). Note 
that values for C° p { 100 ) and C~, 100 . correspond to the tilted 
configuration. 
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The screened charge of the silicon ion induces a coupled 
rehybridization and relaxation: while the silicon neigh- 
bor shell relaxes outward (4.1 Bohr Si-Sixsi distance), the 
carbon second neighbor shell moves inward with a C-Sixsi 
distance of 3.94 Bohr. Apparently, the unfavorable Si-Si 
bonds makes this interstitial site less favorable than the 
SiTC-site. This explains also the observed instability of 
the hexagonal site, which would lie on the line connecting 
adjacent Sixsr and Sixc-sites: The three-fold carbon co- 
ordination cannot counterbalance the energy costs of the 
unfavorable silicon-silicon bond distance until the tetra- 
hedrally carbon-coordinated site is reached. 

In the split-interstitials, as displayed in Fig.|3 the two 
silicon atoms that comprise the interstitial have a mutual 
distance of 4.1 Bohr (Si sp ( 110 )) and 6.1 Bohr (Si sp ( 10 o))- 
The bonds of the Si sp (no)-interstitial are rearranged such 
that each of the silicon atoms has a three-fold carbon 
coordination and a weak cr-bond to its silicon partner. 
Within the band gap we find two non-degenerate deep 
levels below mid-gap (approximately at Ey + 0.4 eV and 
i?v+0.6 eV) that possess anti-bonding character between 
the two silicon atoms. There is also a level close to the 
conduction band. Hence, one may expect that Si sp ^no) 
could exist in the charge states 4 + to 2~. However, ac- 
cording to our calculations only the charge states 2 + , 1 + 
and neutral are realized. In the Si sp ( 10 o) -interstitial both 
silicon atoms are located in the interstitial 3.1 Bohr away 
from the silicon lattice site. This structure corresponds 
to a complex of two Sixc-interstitials and a silicon va- 
cancy. We do not find deep states in the band gap and 
the screening is almost prefect. According to our calcu- 
lations this interstitial is almost unstable. However, it is 
an important intermediate configuration in the migration 
of the Sixc-interstitial. 



B. Carbon interstitials 

A different hierarchy is found for the carbon intersti- 
tials as for the silicon interstitials. In Fig. [S] we have de- 




front side 



FIG. 5: Carbon-interstitials: (a) formation energy for Si-rich 
conditions in comparison to the carbon vacancy (b) geometry 
of investigated interstitial sites. 



picted the stable carbon interstitial configurations. The 
formation energies vs. the Fermi-level are displayed in 
Fig. Eli (for silicon rich conditions) in comparison with 
the formation energy of the carbon vacancy (c.f. Sec. UVfl . 
In relation to Eq. Q we again refer to Si-rich condi- 
tions. For C-rich conditions the formation energy of all 
carbon interstitials drops by 0.61 eV, whereas that of the 
vacancy rises correspondingly. The split-interstitials with 
(lOO)-orientation (C sp (ioo) on the carbon sublattice and 
C sp Si(ioo) 011 the silicon sublattice ) have the lowest for- 
mation energy under all doping conditions. The hexago- 
nal interstitial follows next in the hierarchy. The tetra- 
hedral interstitials and the split-interstitial with (110)- 
orientation are less favorable. For all interstitials deep 
levels exist in the band gap. The outlined hierarchy re- 
flects the preference of carbon to form short bonds in 
the carbon split-interstitial configurations C sp ^ 10 o) and 

C S pSi{100)- 

For C sp (ioo) we find two different configurations de- 
pending on the charge state, as illustrated in Fig. [S] 
In the positive charge states the axis of the carbon 
pair is oriented along the (100) direction. In the neu- 
tral and negative charge states this orientation becomes 
metastable. In the stable configuration the axis of the 
carbon pair is tilted in the corresponding {011}-plane. 
The new stable configuration arises from a pseudo Jahn- 
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Teller distortion of the untilted defect. 

For C^j"/ 1QO j the carbon pair has a (lOO)-orientation. 
The symmetry of the defect is D 2 d- We find a bond dis- 
tance of 2.38 Bohr within the carbon-pair and a silicon- 
carbon bond distance of 3.53 Bohr, which is close to 
the undistorted bond distance in SiC. Along the (100)- 
axis the original tetrahedral arrangement of the silicon 
neighbors is substantially distorted. The bond angle be- 
tween the carbon-silicon bonds amounts to 137.5°. The 
electronic structure is dominated by the sp 2 -like hy- 
bridization of the two carbon atoms. There is a two- 
fold degenerate level (e-representation) within the band 
gap. Besides four carbon-silicon a-bonds, the local- 
ized wave function contains an unhybridized p-orbital at 
each carbon-atom which is oriented perpendicular to the 
plane spanned by the cr-bonds. In the case of Cjp, 100 > 
the partial occupation of the degenerate state invokes a 
Jahn- Teller distortion. The distortion twists the defect 
molecule about its axis by 3.7° so that the silicon neigh- 
bors leave their original positions in the {011}-planes. 
This deformation lowers the energy as it allows a stronger 
interaction of the carbon p-orbital with the silicon or- 
bital. The symmetry becomes D2. Besides the positively 
charged defect also the neutral and negative defect can 
be realized in this twisted configuration. However, in the 
neutral charge state the spin triplet is preferred over the 
singlet by 0.1 eV and the D 2 d-symmetry is retained. 

A tilting of the (lOO)-axis by 31° (in Fig. El in the 
plane Sii-Ci-Si2) transforms the metastable configura- 
tion of C" p ^ 100 ^ and C^ 100 ^ into a stable configuration. 
Alongside with the tilt a simultaneous shift of the atom 
pair occurs, such that the distance of the two carbon 
atoms to S12 become similar (Ci-Si2: 3.7Bohr and C-Si2: 
3.6 Bohr for C° p ^ 100 ^ and C g ~, 100 v). While the distance 
within the pair (2.48 Bohr) only slightly increases, the 
distance Sii-Ci is reduced to 3.2 Bohr. The distances 
Si3-C and Si4-C remain unaffected. The tilt and shift 
introduce a rehybridization of the bonds. Besides the 
carbon-carbon bond, the carbon atom C now binds to 
three silicon neighbors Si2, S13 and S14, while Ci main- 
tains bonds with only two silicon neighbors - a stronger 
one with Sii and a weaker one with Si2- The rebonding 
yields an energy gain of 0.6 eV. The metastable (100)- 
oriented configuration is separated from the tilted ground 
state configuration by a low energy barrier. In the neu- 
tral state the estimated energy barrier is only 0.06 eV. 

The split-interstitial C sp si(ioo) nas a C2 V -symmetry. 
Hence, for this interstitial only non-degenerate levels 
occur within the band gap, which have a similar sp 2 - 
character as C sp (ioo)- There are three levels around mid- 
gap. The lowest one has tr-character. The two higher 
ones have a bonding character with the carbon neigh- 
bors and an admixed p-character on the carbon atom. 
CspSi(ioo) occurs in the charge states 2 + to 2~, with the 
lowest electronic level being always completely occupied 
(c.f. the ionization levels in Tab. QJ. The bond length 
of the atom pair is slightly shorter than the ideal SiC- 




FIG. 6: Geometry of the carbon split-interstitial C a p(ioo) : 
(a) (lOO)-orientation with D2d-/D2-symmetry, (b) tilted 
Csp(ioo)- The C S pSi(ioo) -configuration corresponds to the un- 
tilted C sp (ioo) with the C-Ci-pair replaced by a Si-Ci-pair and 
C-neighbors instead of Si-neighbors. 



bond. This is also true for the bond between the silicon 
atom of the pair and its other carbon neighbors. The 
carbon-carbon bonds with a length of 2.8 Bohr are con- 
siderably shorter. Note that Jahn- Teller distortions are 
absent in this case. However, with the occupation of the 
non-degenerate defect levels the position of the silicon- 
carbon pair shifts along its axis. 



The hexagonal interstitial Ch cx as reflected by its C3 V - 
symmetry has a non-degenerate level above mid-gap and 
a two-fold degenerate level below the conduction band. 
The states originate from a localized p-orbital that is 
oriented along the symmetry axis and from weak sp 2 -like 
bonds with the carbon atoms of the hexagonal ring. In 
the relevant charge states 2 + , 1 + and only the non- 
degenerate p-like level is occupied. As a result the inter- 
stitial carbon is shifted along the hexagonal axis out of 
the center of the ring towards the carbon neighbors (bond 
distance 3.04 Bohr as compared to the distance to the sil- 
icon neighbors of 3.7 Bohr). In the charge state l - the 
occupation of the degenerate level leads to Jahn- Teller 
distortion, that lifts the axial symmetry. As a result, the 
interstitial carbon approaches two of the carbon neigh- 
bors. 



A similar behavior is observed for the tetrahedrally 
carbon-coordinated configuration (Ctc)j which has a 
non-degenerate and a three-fold degenerate level within 
the band gap. Yet, the three-fold coordinated configura- 
tion is energetically preferred over the tetrahedral config- 
uration, since in the latter case the optimization of the 
carbon-carbon bonds results in an unfavorable elonga- 
tion of surrounding carbon-silicon bonds. At the silicon- 
coordinated site favorable carbon-carbon bonds are not 
available. Therefore these configurations as well as the 
(llO)-oriented split configuration are scarcely occupied 
in the thermodynamic equilibrium. 
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IV. VACANCIES AND 
VACANCY- ANTISITE-COMPLEXES 

A. Vacancies 

The properties of the carbon and silicon vacancies (Vc 
and Vsi) in 3C- and 4H-SiC have been discussed in detail 
in the literature pl&SiSSii^ifeiSiil We summarize these re- 
sults as necessary ingredients for the following discussion 
of the vacancy-antisite complexes and the mechanisms of 
the vacancy migration. 

In Figures Eland |31 our results for the formation energy 
of Vsi and Vc are compared to the formation energy of 
the corresponding interstitials for Si-rich conditions. Vc 
is by far more abundant than carbon and silicon intersti- 
tials. This dominance also prevails in stoichiometric and 
C-rich material. Vsi, on the other hand, is less abundant 
than the interstitials in p-type material and dominates 
under intrinsic and n-type conditions. Furthermore, it 
has a higher formation energy than carbon vacancies. 
Hence, the equilibrium concentration of silicon vacancies 
is by several orders of magnitude lower. 

The electronic structure of both unrelaxed vacancies 
derives from the four dangling bonds of the nearest neigh- 
bors. They have T^-symmetry and possess a three-fold 
degenerate one-electron level (^-representation) within 
the band gap. The corresponding non-degenerate a,\- 
level falls below the valence band edge. The occupation 
of the degenerate ta-level either induces a Jahn- Teller 
distortion or leads to the formation of a non-degenerate 
high-spin state. The particular distinction between Vc 
and Vsi is the extended character of the silicon dangling 
bonds of the former versus the strong localization of car- 
bon dangling bonds of the latter^ Driven by the over- 
lap between the silicon dangling bonds, a considerable 
Jahn- Teller distortion of Vc is observedi^iiSi For Vsi, 
in contrast, correlation effects are more importan1*iS44i^ 
than the electron-phonon coupling. As a consequence, 
a multiplet state is formed^^! which leads to the pref- 
erence for a high-spin statesjiSi^ as predicted by DFT- 
calculations^^ for all charge states. 

In Tab.llTlwe have summarized the ionization levels of 
the two vacancies in 3C- and 4H-SiC for the cubic and 
hexagonal site. In the following we judge the realization 
of the different charge state using the experimental band 
gap as noted in Sec. [HJ According to our findings, Vsi 
exists in the charge states 1 + through 2~ in 3C-SiC. In 
4H-SiC also the charge states 3 _ and AT are possible. 
Vgf is not realized as a ground state with the inclusion 
of c harg e state corrections, a finding^, that is verified by 
Ref. l25l Our ionization levels are in good agreement with 
the results of Ref. obtained including charge state cor- 
rections. A slight deviation from the results of Zywietz 
et alml originates from their neglect of these corrections 
and a different definition of the reference for the negative 
charge states. The realization of the charge states 3~ 
and 4~ is also observed in DFT-LSDA calculation in 6H- 
SiC^S using the LMTO-ASA Greens function approach 



as well as the charge state 2 + . In comparison with our 
results in 4H-SiC, the calculated ionization levels agree to 
within 0.3eV including the levels e(3~|2 _ ) and £-(4~|3 _ ). 
We consider this as a good agreement taking the varia- 
tion of the levels among the different polytypes and the 
inequivalent lattice sites into account (it was not stated 
in Ref. |2(] to which site the calculations referred) . 

The carbon vacancy exists only in the charge states 
2 + and neutral in 3C-SiC. An enhanced electron-phonon 
coupling found for V c and V c ~ destabilizes the charge 
states 1 + and 1~. The negative- U amounts to 0.15 eV 
(Vq). However, as the recent identification of Vj by a 
comparison of the experimental hyperfine tensors and the 
calculated hyperfine interaction in 4H-SiG-£ indicates, 
the prediction of the small negative- U maybe artificial. 
In 4H-SiC also the charge states 1~ and 2~ can be real- 
ized apart from the negative- U effect found for Vq . Note 
that the ionization levels of Vc vary between the cubic 
and hexagonal site. This stems from the different ar- 
rangement of the third nearest neighbors and translates 
into a different hybridization of the defect levels. 



B. Vacancy-antisite-complexes and the 
metastability of Vsi 

For the migration two vacancy-antisite complexes are 
important: VsrSic and Vc-Csi, where the antisite is the 
nearest neighbor of the vacancy. These complexes are 
obtained from Vc and Vsi, by a displacement of a sili- 
con/carbon neighbor into the empty lattice site as illus- 
trated in Fig. □ 

We have found that the VsrSic complex does not pos- 
sess a stable configuration in all relevant charge state. 
As a result it decays into the carbon vacancy. The origin 
is a size effect: both, the silicon antisite and the silicon 
vacancy, show a strong outward relaxation^ For the an- 
tisite this relaxation optimizes the silicon-silicon bond 
distance to its neighbors, yet the bond remains strained 
as the ideal value of silicon bulk cannot be achieved. In 
the vacancy-antisite complex the missing silicon neigh- 
bor gives room for further relaxation by moving the an- 
tisite towards the vacancy. Additionally, the energy of 
the dangling bond located at the antisite is reduced by 
interacting with the carbon neighbors of the vacancy. As 
a result the barrier towards the vacant site vanishes. The 
unstable vacancy-antisite complex transforms into a sta- 
ble carbon vacancy. From the above arguments we expect 
that our result for 3C-SiC 2 — is also valid in other poly- 
types. This expectation is supported by a DFT based 
tight-binding calculation for the neutral complex in 4H- 
SiCM 

In contrast to VsrSic the carbon vacancy-antisite 
complex Vc-Csi is stable. In fact, it is even more sta- 
ble than the silicon vacancy in p-type material and for a 
Fermi-level in a mid-gap range. This is shown in Fig. [S] 
Before we discuss the details and consequences of this 
metastability of the silicon vacancy, we briefly turn to 



FIG. 7: Geometry of (a) the carbon vacancy Vc, (b) the silicon vacancy Vsi and (c) the carbon vacancy-antisite-complex 
Vc-Csi- The hop of a carbon neighbor that transforms the silicon vacancy into a carbon vacancy-antisite-complex is indicated 
by an arrow in (b). 




M-p [eV] 

FIG. 8: Formation energy of the silicon vacancy Vsi and the 
carbon vacancy-antisite-complex Vc-Csi for 3C-SiC and 4H- 
SiC. The experimental band edges of 3C- and 4H-SiC are 
indicated by vertical bars. 



TABLE II: Ionization levels of the native vacancies and the 
carbon vacancy-antisite complex in 3C- and 4H-SiC (cubic 
and hexagonal sites are denoted as k and h respectively). 
The results have been obtained using a 216 atom cell (3C) 
and a 128 atom cell including spin-polarization for the silicon 
vacancy. 
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the electronic structure of the complex. 

Similar to the vacancies its localized levels within the 
band gap derive from the four dangling bonds at the 
carbon antisite and the silicon neighbors of the vacancy. 
By the C3 V -symmetry of the ideal defect these dangling 
bonds give rise to two non-degenerate a-levels and a dou- 
bly degenerate e-level. As for the vacancies one a-level 
forms a resonance in the valence band. The other a-level 
lies at Ey+1.8eV and the e-level is found at Ey+2.1eV 
((Vc-Csi) + ). Whereas the character of the localized a- 
levels is dominated by the carbon dangling bond at the 
antisite, the e-level derives from the three silicon dan- 
gling bonds. In the charge states 2 + , 1 + and the a- 
level is occupied. The occupation of the e-level for (Vc- 
Csi) - gives rise to a Jahn- Teller effect, which lowers the 
symmetry to Cih- According to our findings, the charge 
state (Vc-C Si ) 2 - is unstable (c.f. Tab. HJl. In 4H-SiC 
the carbon vacancy-antisite complex shows similar prop- 
erties. For our present purpose it is sufficient to treat the 
complexes aligned along the hexagonal axis. In this case 
the vacancy and the antisite occupy either neighboring 
hexagonal or cubic sites. The larger band gap of 4H-SiC 
-3.3 eV as compared to 2.4 eV in 3C-SiC- enables the 



negative charge states 2 and 3 . 

The metastability of Vsi and its stabilization in n-type 
material occurs in both 3C- and 4H-SiC. As Fig.[H]shows, 
Vc-Csi is more stable than Vsi under p-type conditions 
and for a mid-gap Fermi-level. Its transformation into 
Vc-Csi is associated with a large energy gain (about 4eV 
for a Fermi- level at the valence band). Since the defect 
levels of the Vc-Csi-complex are located above the levels 
of Vsi, the energy gain is reduced by the successive occu- 
pation of defect levels. This leads to a stabilization of Vsi 
over Vc-Csi for n-type conditions. Only a metastability 
was observed^ using special k-points in a 64 atom cell. 
In this case, the formation energy of the silicon vacancy 
is higher in the negative charge states than the value we 
obtain in the 216 atom cell. The silicon vacancy and the 
vacancy-antisite complex have similar formation energies 
for a Fermi-level above mid-gap in 4H-SiC, for the cubic 
complex at /ip=1.8eV and for the hexagonal complex at 
/iF=2.0eV. In 3C-SiC the formation energies match at 
/ip = 1.8eV. The strong interaction between the vacancy 
and the antisite is reflected in the binding energy of the 
complexes. The binding energy ranges between 1.2 eV 
(2+) and l.OcV (1"). 

Using a LMTO-ASA Greens function approach based 
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on DFT and LSDA Lingner et alm& obtain similar re- 
sults for 6H-SiC. A metastability of the silicon vacancy 
is also predicted by their calculations. The Fermi-level at 
which the Vsi and the complex have similar energies is 
located above 2.8 eV. However, in their calculations the 
metastability is more pronounced which leads to a larger 
energy gain of the transformation. According to Fig. 7 
of Ref. |2fJ the energy difference amounts to 4eV (2.6 eV 
for the unrelaxed defects) for the neutral defects. Using 
DFT-based tight-binding calculations the authors obtain 
1.6eV. Our calculations yield 1.1 eV for 3C-SiC, 1.1 and 
1.3 eV for cubic and hexagonal defects in 4H-SiC, which 
is much closer to the tight binding result. For the carbon 
vacancy-antisite complex they obtain similar ionization 
levels for the positive and neutral charge states of the 
unrelaxed complex. Including relaxations a negative-?/ 
effect is found that is not observed in our calculations. 
Yet, as noted above the agreement of the ionization levels 
of the silicon vacancy is good. It indicates that the for- 
mation energy of the silicon vacancy is essentially shifted 
by about 3eV (1.3 eV for the unrelaxed defects) with 
respect to the value of the carbon vacancy-antisite com- 
plex as compared to our calculations. Consequently, the 
Fermi-level at which the silicon vacancy and the complex 
have a comparable formation energy is found closer to 
the conduction band edge as in our case. 



V. MIGRATION 

A. carbon vacancies 

For the migration of the vacancies we have analyzed 
two mechanisms that involve either nearest neighbor hops 
as observed in silicon^ or second neighbor hops, that 
are e.g. relevant for the gallium vacancy in GaAsi^S^fliSi 
The two hops are depicted in Fig. |5| As discussed in 
Sec. I1VI nearest neighbor hops transform the vacancy 
into a vacancy-antisite complex. In successive hops the 
vacancy would proceed leaving behind a chain of an- 
tisites, which is energetically unfavorable. In the ring 
mechanism proposed by Van Vechten^ this is avoided 
by the vacancy passing twice through the same six-fold 
ring. We have shown in Sec. llVl that the silicon vacancy- 
antisite complex is unstable. For the carbon vacancy, 
this finding rules out any migration mechanism based on 
nearest neighbor hops. 

In our analysis of the migration by second neighbor 
hops and the calculation of the migration barriers, we 
follow the procedure in Ref. 36. There it was shown for 
the gallium vacancy in GaAs that the migration mech- 
anism involves the motion of the hopping gallium atom 
(a C-atom in Fig. 0) and its common neighbor with the 
vacancy (the Si-neighbor in Fig.^Jl. At the saddle point, 
the hopping atom passes through the plane spanned by 
the four atoms (Ci,. . . in Fig. |3J) that we refer to as 
the gate below. In order to accommodate the strain of 
this configuration it is favorable for the nearest neigh- 




FIG. 9: Migration of the carbon vacancy Vc: (1) initial 
nearest neighbor hop, (2) migration on the carbon sublattice 
by second neighbor hops. 

TABLE III: Vacancy migration: the migration barrier and 
geometry at the transition state for the migration by second 
neighbor hops and the transformation of Vsi into Vc-Csi (ini- 
tial nearest neighbor hop). See the text for details. 
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bor to relax towards either of the vacant lattice sites. 
Thus two saddle points and correspondingly two migra- 
tion paths exist for the hop - one, where the neighbor of 
the vacancy passes through the gate prior to the second 
neighbor and one where the latter passes first. Following 
Ref. I36L we introduce two reaction coordinates £c and 
£si for the hopping carbon atom and the common silicon 
neighbor by 

£c/Si = ^R-c/Si _ 4 5^ ^ C i\ ' eno ' ( 7 ) 

where R-c/si f°r £c and £sb respectively, refers to the 
coordinates of the carbon or silicon atom and the Rq ; 
are the coordinates of the four carbon atoms comprising 
the gate; eno is the unit vector of the direction connect- 
ing the initial and final lattice site of the vacancy (here 
the (llO)-direction is chosen). Instead of evaluating the 
full potential energy surface for these two coordinates, 
we consider only the region in the vicinity of the transi- 
tion state. For a given £c the total energy is minimized 
with respect to £si and the other atomic coordinates. We 
find indeed the behavior we have discussed above. The 
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transition state lies by 0.2 eV lower than the configura- 
tion where the carbon atom and its silicon neighbor pass 
through the gate at the same time. The barriers are ob- 
tained using the 216 atom cell together with the T-point, 
which also allowed for the proper inclusion of Jahn- Teller 
distortions in the initial configurations. The results for 
the migration barrier and the geometry at the transi- 
tion state are given in Tab. IIIII Using a 64 atom cell 
and special k-points (2x2x2 Monkhorst-Pack mesh), 
the calculated barriers agree within 0.3 eV. We find rela- 
tively high migration barriers between 5 eV (Vq + ), 4.1 eV 
(Vq) and 3.5eV (Vq). The migration barrier in n-type 
material and with a Fermi-level around mid-gap (Vq) is 
much smaller than in p-type material (V^ + ). This stems 
from the different occupation of bonding defect levels at 
the transition state. Besides the bond between the com- 
mon Si-neighbor and Sio (c.f. Fig. 01, the remaining three 
dangling bonds at the opposite vacant site contribute to 
the defect levels in the band gap. With the occupation of 
the lowest level the corresponding neighbor (Sii) relaxes 
towards the vacant site, thus reducing the distance to the 
silicon neighbor and the migration barrier. 



B. silicon vacancies 

In case of the silicon vacancy the migration by second 
neighbor hops competes with the transformation Vsi <-> 
Vc-Csi- Particularly in p-type material, the metastabil- 
ity of the silicon vacancy will have consequences for its 
migration. However, a migration of the silicon vacancy 
solely based on nearest neighbor hops is impossible. This 
is due to our finding that the second configuration in such 
a mechanism, the Vsi-Sic-Csi-complex, is unstable in all 
relevant charge states. In order to understand the conse- 
quences of the metastability, the barriers for the second 
neighbor hop and the transformation Vs; <-> Vc-Csi have 
to be compared. 

The analysis of the second neighbor hop proceeds as 
detailed in the previous section. This mechanism has es- 
sentially the same features as in the case of the carbon 
vacancy. Now the second silicon neighbor and the near- 
est carbon neighbor take the role of the respective carbon 
and silicon atom in the discussion of Vc- The calcula- 
tions have been performed within the 216-atom cell us- 
ing the T-point and including spin polarization. For the 
64-atom cell and a special k-point sampling we obtain 
within 0.2 eV the same barriers. Our results for the bar- 
riers and the geometry at the transition state are listed in 
Tab. 11111 At the transition state the displacement of the 
carbon neighbor out off the gate is less pronounced as in 
the case of the carbon vacancy. For the negative vacancy 
it almost vanishes. We explain this by the rather small 
extend of the carbon orbitals and the resulting weak in- 
teraction of the carbon neighbor with the other carbon 
dangling bonds. The migration barrier for this mecha- 
nism varies between 3.6 eV (V^) and 2.4 eV (V|r). 

Our analysis of the nearest neighbor hop parallels the 



approach to the second neighbor hop. We introduce the 
reaction coordinate £nn which measures the distance of 
the hopping carbon atom from the center of mass of its 
three silicon neighbors. It is given by 

£nn = ( Rc - g ^ Rsij J • em , 

where Rc is the coordinate of the hopping carbon atom 
and Rgii are the coordinates of its three silicon neighbors; 
em is the unit vector in the direction of the hop (here 
chosen as the (lll)-direction). Note that this approach 
allows for a migration path that deviates from the (11 In- 
direction. 

As detailed in Sec. II VI Vsi and Vc-Csi are found in dif- 
ferent charge states for a given Fermi-level. We assume 
that the charge state of the defect is preserved during the 
transformation. The necessary electron transfer to obtain 
the equilibrium charge state occurs in the final configura- 
tion. Harrison has shownS for the case of the interstitial 
in silicon that this assumption is reasonable. This ap- 
plies also to the transformation of V S J~ — > (Vc-Csi) 2- 
and (Vc-Csi) 2+ ^V 2 .^, where the final configurations 
(Vc-Csi)~ and Vgj are the ground state configurations. 

The energy barriers and the values of the reaction co- 
ordinate £nn at the transition state are listed inTab.Hnl 
The barrier for the transformation Vsi — > Vc-Csi mod- 
erately depends on the charge state. It is by 0.8 eV lower 
for the positive vacancy than for the doubly negative 
charge state. The reverse transformation Vc-Csi — » Vsi 
shows a much stronger charge state dependence. A bar- 
rier of 6.1 eV is found for (Vc-Csi) 2+ , that lowers down 
to 2.4 eV for (Vc-Csi) - - During the transformation 
the hopping carbon atom moves away from the (1 In- 
direction. In this way it optimizes the number of bonds 
maintained with its silicon and carbon neighbors. A 
path along the (11 Indirection has a much larger energy 
barrier, due to a simultaneous breaking of three bonds 
with the silicon neighbors. Only this latter path was 
obtained^ using the 64 atom cell with a special k-point 
sampling and was reproduced using the 216 atom cell. 
However, this path is not important. 

For all charge states the transformation Vsi ^Vc-Csi 
mainly occurs in a low-spin state of the vacancy. Typ- 
ically much higher barriers are found for the high-spin 
state, except for Vg ; where the high-spin state has a 
slightly lower barrier (given in braces in Tab. 1111(1 . Note 
that the ground state of the vacancy-antisite complex is 
always a low-spin state for all charge states. Thus the 
reverse transformation Vc-Csi — >Vsi should occur into 
the low-spin state. 

In equilibrium (with a given Fermi- level), Vsi and Vc- 
Csi are in different charge states. Provided the temper- 
ature is sufficiently low and the transformation is a rare 
event, the charge state can equilibrate at the final po- 
sition before the reverse process takes place. Only at 
extreme temperatures we expect this picture to break 
down. Consequently, energy barriers for different charge 
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states should apply for the two processes. In p-type 
and compensated material the barrier for the transfor- 
mation V Si ^Vc-Csi is with 1.9 eV (V s +) to 2.4 eV 
(Vgj) much lower than the barrier of the reverse trans- 
formation of 6.1 eV ((V c -C Si ) 2+ , mf < 1.24 eV) to 3.5eV 
((Vc-Csi)°,^F < 2.19 eV). In n-type material (i.e. for 
> 2.19) this situation is reversed, when (Vc-Csi) - 
becomes the ground state. Now the transformation of 
(Vc-Csi)~ into has a higher probability than the op- 
posite process Vgj~ — >(Vc-Csi) 2_ . 

A comparison between the migration barriers of the 
silicon vacancy and the transformation barriers should 
tell whether the silicon vacancy is trapped as a Vc-Csi- 
complex or whether it freely migrates. In n-type material 
Vgj" is stable (i.e. for > 1.76eV). The migration bar- 
rier (2.4 eV) is slightly smaller than the transformation 
barrier (2.7 eV) to the metastable Vc-Cgi-complex. For 
Hf > 2.19 eV (Vc-Csi)~ returns by a reverse transforma- 
tion (2.4 eV) via Vg; to V|f with an even higher proba- 
bility. For fi F between 1.79eV and 2.19eV (V c -Cgi)° is 
relevant. Even though the barrier of the transformation 
Vc-Csi — >Vsi (3.5 eV) is larger, the migration of Vs; is es- 
sentially not hindered as Vg~ is the ground state. Thus, 
in n-type material, the vacancy should freely migrate. 

On the other hand, for a Fermi- level below 1.7 eV the 
Vc-Csi-complex is the ground state. Before a migration 
of Vs; may take place, it has to be preceeded by a trans- 
formation Vc-Csi — *Vsi- This gives rise to an additional 
barrier. For a Fermi-level below mid-gap this barrier be- 
comes large, i.e. 4.1 eV for (Vc-Csi) + and 6.1 eV for 
(Vc-Csi) 2+ ■ This additional barrier allows a migration of 
Vc-Csi via Vsi only at elevated temperatures, when the 
transformation Vp-Csi — +Vsi is thermally activated. At 
low temperatures this migration is kinetically hindered. 

Instead, another migration channel is opened: the Vc- 
Csi complex can dissociate at the expense of the binding 
energy. However, using the migration barrier for Vc be- 
tween 3.5 eV (Vp) to 5eV (Vq + ) plus the binding energy 
as an estimate for the dissociation barrier, we conclude 
that the dissociation and the reverse transformation have 
similar probabilities. In p-type and compensated mate- 
rial {up < 1.71 eV) the migration of silicon vacancies is 
suppressed for moderate temperatures. 



C. carbon-interstitials 

Apparently, the energetically lowest sites are the most 
relevant ones for the migration of the carbon intersti- 
tials. As discussed in Sec. IIIII these are the two split- 
interstitials C sp (ioo) and C sp si(ioo)- The migration of the 
energetically more favorable C sp ^ 1 oo)-interstitial may pro- 
ceed by two alternative routes: (a) by nearest neighbor 
hops between adjacent carbon and silicon lattice sites in 
the sequence C sp(100 ) -*C spSi (ioo) ^C sp(100) and (b) by 
second neighbor hops of the interstitial between neigh- 
boring carbon lattice sites. 

We illustrate the migration event by nearest neigh- 



TABLE IV: Migration of carbon interstitials: Migration bar- 
riers E m in eV for nearest neighbor hops C sp (ioo) <_i 'C a pSi(ioo> 
and second neighbor hops. The migration of Cg p ( 100 ^ and 
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bor hops in Fig. ^3 Taking C sp ( 10 o) (Fig. ITUk ) as the 
initial configuration, the final configuration is C sp gj{ 100 ) 
(Fig. 110b). which is located at the nearest neighbor sil- 
icon site. The saddle point search yields the transition 
state shown in Fig. HUb . The migration proceeds by a 
concerted motion of the dumbbell, such that one carbon 
atom approaches the neighboring silicon atom while the 
other returns into the position on the lattice. At the same 
time, the silicon neighbor moves out of this site to give 
room for a new carbon-silicon dumbbell. We have con- 
ducted a saddle point search in the charge states 2 + , 1 + , 
neutral and 1~ using a 64 atom cell and a special k-point 
set. For the hop from C sp (ioo) to C sp si(ioo> we nn d a bar- 
rier of 1.7 eV and 0.9eV for the charge states 2 + and 1 + . 
The barriers of the second hop from C sp si(ioo) to C sp (ioo) 
amount to 0.7 eV and 0.2 eV respectively and are lower 
than those of the first hop. In the neutral and negative 
charge state the hop starts from the tilted configuration. 
The path is as discussed above. The barriers amount to 



p(100) 



Similar barriers are found 



0.5 eV and 0.2 eV for C, 
in the case of CV ^ 100 ^. 

The migration of C sp (ioo) by second neighbor hops can 
proceed in several different ways. Each of the carbon 
atoms of the dumbbell is surrounded by six neighboring 
carbon lattice sites. While one atom migrates to one of 
the neighboring carbon sites, the other moves back onto 
the lattice. We have investigated the hop of the (001)- 
oriented split-interstitial in the (llO)-direction. The mi- 
gration proceeds in a similar manner as the nearest neigh- 
bor hop, except that now the destination is the adjacent 
carbon lattice site. The migration barriers we obtained 
with the saddle point search are listed in Tab. IIVI They 
are comparable to the nearest neighbor hop. Only for 
Cwioo) we find a lower barrier for the second nearest 
neighbor hop (1.4 eV) than for the nearest neighbor hop 
(1.7 eV) 

The relatively small migration barrier contrasts the 
much larger energy difference found between the split- 
interstitials and the tetrahedrally coordinated sites. As 
at these sites carbon bonds are 
Consequently, a rele- 



discussed in Sec. IIII Bl 
established at high energy costs, 
vant migration of carbon interstitials via these sites is 
unlikely. 

The hexagonal interstitial is 1.4 eV (CgL.) and 1.6 eV 
( C Hex) higher in energy than C S p Si (ioo). This energy dif- 
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ference is only slightly lower than the migration barrier 

for C s^Si<100) and eVen hi S hcr in tllC CaSe ° f C spSi(100>- 

Therefore we expect a higher energy barrier for a migra- 
tion via the hexagonal sites, particularly as the mecha- 
nism should initially proceed similar to the migration via 
split-interstitial sites to avoid extensive bond breaking. 



D. silicon-interstitials 

In Section ITfTAl we have shown that the two most fa- 
vorable interstitial configurations are the tetrahedrally 
carbon-coordinated site and the (llO)-oriented split- 
interstitial. The first configuration is favored in p-type 
material and the latter is clearly dominating in n-type 
material. The migration between the energetically lowest 
sites may involve the other configuration as an interme- 
diate site or may proceed by direct hops. As discussed 
previously, we assume that the charge state is not altered 
during a hop. Only in the intermediate configurations a 
charge transfer is allowed. 

In p-type material all migration paths begin and end 
at the SiTC site. Two different types of migration 
mechanisms are discussed here and depicted in Fig. 1771 
(a) a kick-out mechanism that proceeds via the split- 
interstitials Si sp (100) and Si sp mo) as intermediate sites 
and (b) a direct migration between the tetrahedral in- 
terstitial sites Sixc and Sixsi- The kick-out mechanism 
consists of two steps. In the first step the tetrahedral 
interstitial moves towards a silicon neighbor and kicks it 
out of its lattice site. A split-interstitial is formed as an 
intermediate site. Next, the second atom proceeds to- 
wards the closest tetrahedral interstitial site. We have 
determined the migration paths and the energy barriers 
for the two processes by a saddle point search using a 64- 
atom cell and special k-points (2x2x2 Monkhorst-Pack 
mesh). While the path via the Si sp (ioo) -configuration 
proceeds entirely along the (lOO)-direction, the path via 
the Si sp ^no) interstitial is curved. The transition state 
in both cases lies close to the split-interstitial configura- 
tion. The energy barriers for the kick-out and the kick-in 
process are listed in Tab.lVl The kick-out barrier via the 
Si S p (ioo) interstitials amounts to 3.5 eV and is only 0.1 eV 
higher than the barrier obtained via the Si sp (no) intersti- 
tial. The results show that the two intermediate intersti- 
tial configurations are not very stable as the barriers for 
the second event amount to 0.03 eV and 0.05 eV, respec- 
tively. We expect that the interstitial passes through the 
intermediate site without any significant delay. There- 
fore a charge transfer that could convert the Si sp { 110 \- 
interstitial from the transient 4 + -state into 2 + or 1 + (the 
ground state in p-type material) seems unlikely. 

The migration between the tetrahedral interstitial sites 
SiTC and Sixsi proceeds directly through the interstice. 
The migration path passes through the hexagonal con- 
figurations along the axis connecting the initial and fi- 
nal configuration. The transition state lies close to the 
ideal position of the hexagonal site. The saddle point 



TABLE V: Migration of silicon interstitials: Migration bar- 
riers Em in eV for different migration paths. The migration 
paths are shown in Fig. 1111 The charge state is 4 + unless 
indicated otherwise. See text for details. 
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search yields an energy barrier of 3.5 eV for the hop 
Sixc — ^Sixsi- For the hop Sixsi ~~ >Sixc w ^ obtain a 
barrier of 0.3 eV. 

Since for Sixsi the charge state 4 + is not the ground 
state, a charge transfer can change it with some proba- 
bility. Then the interstitial could alternatively migrate 
via the split-interstitial Sr^~/ 110 » with a low barrier (c.f. 
Tab. IV|) . However, this route is not very likely Instead 
the Si T gj (or Si T g;) would move to the Sixc _ site with an 
even lower barrier. 

Incompensated and n-type material, we consider the 
kick-out migration via the split-interstitials Sig p / 110 v. As 
one of the silicon atoms of Si sp ( 110 ) hops towards the 
next silicon site it kicks out the silicon atom at this site 
in the corresponding (llO)-direction to form a new split- 
interstitial. For this process we find a migration barrier 
of 1.4 eV. For the positive split-interstitial we obtain a 
barrier of 1.0 eV. The orientation of the split-interstitial 
can change by a rotation around the (100)-axis. For ex- 
ample, the rotation around the (OOl)-axis changes the 
orientation from (110) to (110). It is associated with an 
energy barrier of 1.3eV (Sig p ( 110 ^). 

In order to understand whether the kick-out migration 
via Sig p ( 110 ^ is the dominating mechanism, we have to 
compare its activation energy with that of other possible 
mechanisms at intrinsic and n-type doping conditions. In 
comparison to a migration via the interstitial Si T g j7 we 
find that this configuration for /ip above 1.3 eV already 
lies higher in energy than the saddle point of Si sp ( 110 ). 
Furthermore, according to our estimates also Sixc and 
Si S p (ioo) should be correspondingly higher in energy for 
a Fermi-level close to the conduction band. Thus we 
conclude that the dominating mechanism under intrinsic 
and n-type conditions is the migration via the Si sp ( 110 )- 
interstitial. 



14 






Si S p<110> 



FIG. 11: Migration of silicon interstitials: (a) kick-out migration of SiTC via the split-interstitial Si sp (iio>, (b) migration via 
the tetrahedral interstitial sites Si-rc and Slrsi, and (c) jumps between adjacent Si sp /uo) configurations. 



VI. DISCUSSION 

Our analysis of the mobile intrinsic point defects in 
the Sections IIIII and llVl and their migration mechanisms 
in Section Ivl reveals a complex dependence of the abun- 
dance and the migration on the doping conditions. In the 
following we discuss the implications of our findings for 
the self- and dopant-diffusion. The mechanisms that en- 
ter the discussion are summarized schematically in Fig.Q] 
and the associated energy barriers are listed in Tabs. lTTll 
IIVI and This microscopic picture of the point defect 
migration has also implications for the annealing mecha- 
nism of vacancies. This aspect is treated, together with 
an analysis of Frenkel pairs, in a separate articled 

A. Self diffusion 

The self diffusion constants of carbon and silicon, Dc 
and Z?si, are given by the sum of different diffusion mech- 
anisms, in particular of mechanisms based on vacancy 
and interstitial migration. For instance, the carbon self 
diffusion constant is given by 

/; <- J2 D 7 a • • ••• ( g ) 

i J 

where Df c and D\' are the diffusion constants of va- 
... 

cancy and interstitial mechanisms, respectively, and i, j 
denote the different migrations paths of the mobile defect 



in a specific charge state. Each contribution is propor- 
tional to the concentration of the mobile defect Cj (c.f. 
Eq. (JTJ) and its diffusivity di, as given by the formation 
energy E\ and the migration energy barrier Ef 1 , respec- 
tively, 

D, = A°e,p p (A *;g + £ " ) , (9) 

where D® is a constant prefactor that is determined by 
geometrical and entropy factors. The activation energy 
Qi = El(Afi, /xp) + Ef 1 that determines the exponential 
factor essentially characterizes the contribution of each 
migration mechanism. The factor may vary over several 
orders of magnitude for different mechanisms, whereas 
the prefactors D® typically have a similar order of mag- 
nitude. Thus a comparison of the activation energies of 
different diffusion mechanisms reveals the dominant con- 
tribution to self diffusion for different doping and growth 
conditions that enter Q via Afi and fip. In Table IVTI the 
activation energies for vacancy- and interstitial-mediated 
self diffusion are listed. Q is given as an explicit function 
of /ip and A^i. At Si-rich conditions A/z is and assumes 
a value of —Hf lC (— 0.61eV in the present calculations) 
at C-rich conditions. For the charged defects the range of 
is given in which the particular charge state is stable. 
The variation of Q with the Fermi-level can be substan- 
tial, i.e. for Si^Q it amounts up to 5eV. For a typical 
diffusion temperature of 1000 K the change of AQ ~ 1 eV 
translates into a variation of the diffusion constant over 
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TABLE VI: Activation energy for the self diffusion as mediated by vacancies and interstitials (cf. text for details). 



Defect Q(fi c ,An) [eV] ^ F [eV] 

explicit Si-rich C-rich range 

Vq + 6.6 eV +2/ZF -A/x 6.6-9.0 7.2-9.6 0.0-1.2 

V° c 7.3 eV -An 7.3 7.9 1.2-E C 

C sp(ioo> 6 - 7eV + 2 ^ F + A / i 6.7-7.8 6.1-7.2 0.0-0.6 

C sp<ioo> 6 - 6eV +A'f 7.2-7.9 6.6-7.3 0.6-0.8 

Csp(ioo) 7.7 eV +A/i 7.2 6.6 0.8-1.8 

C- (100> 9.1 eV -a*f +Afi 7.3-6.7 6.7-6.1 1.8-E C 

Exp." 3C 8.72eV±0.14eV undoped 

Exp. 6 4H 7.14 eV± 0.05 eV intrinsic 

Exp. 6 4H 8.20 eV± 0.08 eV n-type 



a Ref. 1^3 experiments on poly-crystalline 3C-SiC samples. 
6 Ref. IgflOi experiments on 4H-SiC single crystals. 



Defect Q(fj. u Afi) [eV] /i F [eV] 

explicit Si-rich C-rich range 

Vgf 12.2 eV— fiF +Afi 11.6-10.4 11.0-9.8 0.6-1.8 

Vgr 13.1 eV -2 fi F + A fi 9.6-7.7 9.0-7.1 1.8-E C 

S4c 7.8eV+4^ F -A^ 7.8-12.6 8.4-13.2 0.0-1.2 

7.7eV+4/i F -A^i 7.7-12.5 8.3-13.1 0.0-1.2 

Si° p(110) 9.9 eV -Afi 9.9 10.5 1.1-2.4 

Exp." 3C 9.45 eV ± 0.05 eV undoped 

Exp. 6 4H 7.22 eV± 0.07 eV intrinsic 

Exp. 6 4H 8.18eV±0.10eV n-type 



five orders of magnitude. This striking dependence on the 
doping conditions is known as the Fermi- level effect^ 
Table IVII also includes the activation energies deduced 
from recent self diffusion experiments for polycrystalline 
3C-SiC 7 - s and for 4H-SiC single crystalsAifi To the best 
of our knowledge, these are the first experiments that 
quantitatively analyzed the self diffusion constant. A 
radiotracer diffusion technique was used to obtain the 
diffusion profiles. 

For the carbon self diffusion our results show a simi- 
lar importance of the vacancy-mediated and interstitial- 
based self diffusion in Si-rich material. There is no clear 
dominance of either mechanism. In C-rich material, on 
the other hand, the split-interstitial mediated self diffu- 
sion dominates for all doping conditions. Since in com- 
pensated and n-type material the split-interstitial and 
the vacancy are neutral, the Fermi-level effect is absent. 
Only in p-type material a variation of the activation en- 
ergy about 1 eV is present. 

The silicon self diffusion in p-type and compensated 
material is mainly dominated by the silicon interstitials. 
In p-type material the migration is mediated by the tetra- 
hedral carbon-coordinated interstitial Sixc j both, under 
Si-rich and C-rich conditions. Here the metastability of 
the silicon vacancy suppresses its participation in the 
self diffusion. In Si-rich compensated material the sili- 
con spit-interstitial Si sp ^ 110 ) is relevant. Only in n-type, 
Si-rich material a dominant contribution of the silicon 
vacancy Vsi should be observed. The larger abundance 
of silicon vacancies under C-rich conditions leads to its 
dominance already in compensated material. The varia- 
tion of the activation energy of the silicon self diffusion 
is more pronounced than for the carbon self diffusion. A 
Fermi-level effect is predicted under p-type and n-type 
conditions. The variations amount to 2eV, both, in C- 
rich and Si-rich material. 

Experimentally the silicon and carbon self diffusion 



was investigated in nominally undoped poly-crystalline 
3C-SiCA& For the carbon self diffusion an activation 
energy of 8.72 eV was obtained. The observed silicon 
self diffusion was described 8 by an activation energy of 
9.45 eV. From the description of the experiments it is not 
obvious whether they were carried out in C-rich or Si-rich 
conditions. Therefore we compare the experimental re- 
sults with our values for C-rich and Si-rich conditions 
for an intrinsic Fermi-level. For the carbon self diffusion 
we obtain an activation energy of 7.3 eV (Vq, Si-rich) 
and 6.7 eV (C sp (no), C-rich), respectively. The experi- 
mental value exceeds our theoretical results by 1.4 eV- 
2.0 eV. This discrepancy cannot be explained in terms 
of the missing information about A/i or a variation of 
fip. Our results indicate a value of at most 8eV. Only if 
the interstitial diffusion was completely suppressed dur- 
ing the experiment and the sample was rather p-type 
than intrinsic, the activation energy could be explained 
by that of . Here the poly-crystalline nature of the 
material may also affect the defect equilibrium, i.e. the 
grain boundaries may act as sinks for interstitials. For 
the silicon self diffusion we obtain an activation energy 
of 9.9eV (Si sp(110> , Si-rich) and 10.4eV (Vgj, C-rich) in 
intrinsic material. In this case the experimental value is 
in good agreement with our results. 

Even though our investigation of the defect migration 
was conducted for 3C-SiC, we expect that our results 
may provide a first insight into the diffusion processes in 
other polytypes. This expectation is supported by the 
findingi2iS4i2^ of comparable energetics and bonding of 
vacancies in 3C-SiC and 4H-SiC. For example, similar 
transformation barriers are found for the neutral silicon 
vacancy in 3C-SiC (this work) and 4H-SiCi2i The main 
difference is the larger energy gap in 4H-SiC. Intrinsic 
conditions, thus, apply in 4H-SiC to fip ~ 1.7 eV. A com- 
parison of our results with the activation energy of the 
carbon self diffusion in 4H-SiC single-crystals shows a 
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much better agreement than for the poly-crystalline 3C- 
SiC. However, for the silicon self diffusion our prediction 
of about 9 eV substantially deviates from the experimen- 
tal result of 7.2 eV, unless we place the Fermi-level above 
mid-gap at around 2.3 eV. rt-type 4H-SiC is not covered 
by our results, since the relevant charge states of the 
vacancies and interstitials in 4H-SiC are not present in 
3C-S1C. 

In the light of our expectation the large variation of 
the experimental activation energies between 3C-SiC and 
4H-SiC is surprising. Besides the poly-crystalline nature 
of the cubic material, an incomplete point defect equi- 
librium or the preferential injection of either vacancies 
or interstitials at the interface of the tracer material and 
SiC may be responsible. Here we note that the investi- 
gation of the aluminum and boron diffusion^ indicate 
a slower silicon self diffusion than observed in the tracer 
diffusion experiments. This finding inevitably affects the 
given prefactor and activation energies. 



B. Dopant migration 

In SiC, dopant atoms may substitute for carbon or sil- 
icon. For example, aluminum and phosphorus are known 
to preferentially bind to the silicon sublattica^ and ni- 
trogen prefers carbon sites^ Boron is predicted to bind 
either to the carbon or silicon sublattic o 56 ' 57 depending 
on the stoichiometry of the sample (C/Si-ratio). The mi- 
gration of these dopants is mediated by mobile intrinsic 
defects. A mobile vacancy and a substitutional dopant 
atom may form a mobile dopant-vacancy-complex. An 
attractive interaction within the mobile pair enables mi- 
gration until the pair eventually dissociates. Mobile in- 
terstitials may initiate a migration of the dopant by a 
kick-out reaction. 

In equilibrium, the dominating mechanism of dopant 
diffusion depends not only on the involved migration bar- 
riers but also on the abundance of the involved intrinsic 
defects. A transient enhanced dopant diffusion that was 
recently observed for boroni£ and aluminum* is triggered 
by an excess concentration of intrinsic defects, generated 
for example by the dopant implantation. Assuming that 
a similar amount of interstitials and vacancies is gen- 
erated, the migration barriers for different mechanisms 
are the relevant quantities to compare. Naturally, for 
dopants on the carbon and the silicon sublattice different 
migration mechanisms apply. Even though the migration 
mechanism should certainly be dopant-dependent, some 
general conclusions about the dopant diffusion can be 
drawn from our results for the diffusion of vacancies and 
interstitials. 

First of all, the availability of the mobile defects de- 
pends on the doping conditions. The discussion in Sec. II VI 
showed that a mechanism based on silicon vacancies is 
not available in p-type material. The metastability of Vsi 
in p-type material should affect dopant complexes with 
a silicon vacancy in a similar way. Thus only the migra- 



tion of a donor complex that involves the carbon vacancy 
should be relevant. In such a complex the dopant and 
the vacancy may be nearest neighbors or second nearest 
neighbors, when the dopant substitutes for silicon or car- 
bon, respectively. Whereas the migration of the second 
neighbor pair proceeds entirely on the carbon sublattice 
by second neighbor hops, the migration of the nearest- 
neighbor pair involves an initial exchange of the dopant 
and the vacancy and further migration of the vacancy 
around the dopant i£ s This may only be accomplished by 
second neighbor hops. A mechanism based on nearest 
neighbor hops fails as the necessary Vsi-Sic-complexes 
are unstable. In n-type material a mechanism based on 
silicon vacancies becomes available. With the same ar- 
guments as above a migration will involve only second 
neighbor hops of the vacancy. 

Consequently, the diffusion of acceptors on the sili- 
con sublattice most likely is mediated by silicon inter- 
stitials in p-type material. More general conclusions on 
the role of the mechanisms based on carbon and silicon 
interstitials are not possible. Here we only note that 
carbon-interstitials and silicon-interstitials (in compen- 
sated and n-type material) have lower migration barriers 
than the corresponding vacancies. The interstitial medi- 
ated dopant migration may be limited by the availability 
of interstitials and the barriers associated with the mi- 
gration of the dopant interstitial. 

These arguments are supported by our findings for 
boron££ and by recent experiments on the diffusion of 
boronli2i2£ and aluminum. 4 Our calculations reveal that 
a migration via the nearest-neighbor complex Bgi-Vc is 
unlikely The barrier of the initial transformation has 
a similar magnitude as the transformation of Vc-Csi 
into Vsi in p-type material. Furthermore, the nearest- 
neighbor complex Bsi-Vc is practically unstable in p- 
type material. The migration of the second-neighbor 
complex Bc-Vc consists of three hops and involves mi- 
gration barriers in p-type material of as much as 5.6 eV, 
which exceeds the barrier for the carbon vacancy. We 
find much lower migration barriers for a kick-out mech- 
anism by silicon interstitials (Bsi) or carbon intersti- 
tials (Be) and a subsequent migration of the boron- 
interstitials. Though earlier experiments 6 favored a va- 
cancy related mechanism, the recent quantitative mod- 
eling of boron diffusion experiments^^ attribute the 
diffusion to a kick-out mechanism based on silicon in- 
terstitials. A similar conclusion was suggested by co- 
implantation experiments^ Here the transient enhanced 
diffusion observed after boron implantation was sup- 
pressed by carbon co- implantation and elevated by silicon 
co-implantation. 



VII. SUMMARY AND CONCLUSION 

A picture of the defect migration has been derived from 
theoretical investigations based on an ab initio method 
within the framework of DFT. The analysis of the mi- 
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croscopic structure, the abundance and the migration 
mechanisms of interstitials and vacancies contribute to 
this picture. A strong influence of the doping condi- 
tions is reported for the mobile silicon defects. Among 
the silicon interstitials two different interstitial configura- 
tions are relevant, the tetrahedrally carbon-coordinated 
interstitial in p-type material and the split-interstitial in 
compensated and n-type material. This affects the in- 
terstitial migration in p-type material: the migration of 
the tetrahedrally carbon-coordinated interstitial is slower 
than the migration of the split-interstitial. A metasta- 
bility of the silicon vacancy occurs in p-type and com- 
pensated material. It transforms into the more stable 
carbon vacancy-antisite complex. This transformation of 
the silicon vacancy is more likely than its migration. The 
thermal activation of the reverse process determines the 
mobility of the vacancy-antisite complex. A central re- 
sult is the finding that interstitials are more mobile than 
vacancies. 

The importance of interstitials and vacancies in the self 



and dopant diffusion is determined by their abundance 
and diffusivity. The discussion of the two factors showed 
that the interstitials play a significant role in the self 
diffusion. For instance, silicon interstitials dominate the 
silicon self diffusion in p-type material. The dominance of 
carbon interstitials in the carbon self diffusion stems from 
an over-compensation of their relatively low abundance 
by their higher mobility. For the dopant diffusion, similar 
qualitative conclusions are drawn and outlined for the 
case of boron, for which theoretical calculations indicate 
an important role of the interstitial-mediated diffusion. 
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